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Abstract. We develop an analytical model for describing the dynamics of a donor-based 
charge quantum bit (qubit). As a result, the quantum decoherence of the qubit is analytically 
obtained and shown to reveal non-Markovian features: The decoherence rate varies with 
time and even attains negative values, generating a non-exponential decay of the electronic 
coherence and a later recoherence. The resulting coherence time is inversely proportional to 
the temperature, thus leading to low decoherence below a material dependent characteristic 
temperature. 
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1. Introduction 

Implementations of quantum bits (qubits) in solid-state systems represent some of the most 
promising candidates for quantum-computation devices. The constituent two-level systems 
have been experimentally realized in a variety of systems, including quantum dots [1], 
impurities in diamond [2], and superconductors [3] among others. The implementation of 
solid-state quantum-information devices has been object of increasing interest in the last 
decade [4-6]. 

Like other types of implementations, solid-state systems also suffer from the inevitable 
interaction with their environment, which leads to a loss of coherence of individual qubits 
and of entanglement between them. Whereas the decoherence of many physical systems 
can be well approximated by a Markovian process, a few prominent examples exist where 
the coupling of the system with its environment reveals a memory effect, leading to non- 
Markovian behavior and the possibility of recoherence at later times. Among these are spin 
echo [ ] — which has been recently observed with phosphorus donors in silicon [8], photon 
echo [9, 10], and collapse [11] and revival [12-14] in the spin-boson Jaynes-Cummings 
interaction [15]. Collapse and revival may occur in various physical systems, such as with 
atoms interacting with cavity electromagnetic radiation [16], in the vibronic dynamics of 
electromagnetically trapped ions [ 1 7], or in molecular vibration [18, 19]. However, also other 
systems may show a non-Markovian dynamics in the decay of their coherence. In particular, 
numerical calculations have shown that donor-based charge qubits interacting with phonons 
reveal a non-Markovian decay of electronic coherence [20]. 

Donor-based charge qubits consist of a single electron shared by two impurities in 
the semiconductor host. Besides being potentially scalable, their manufacturing is largely 
facilitated by the technological base provided by the existing semiconductor industry. One 
of the main sources of decoherence in this system is the coupling of the electron to acoustic 
phonons. This decoherence mechanism has been studied numerically for mK temperatures 
[20] and analytically for T = OK [21]. In the former reference, a non-exponential decay of 
coherence has been numerically obtained, which indicates a non-Markovian behavior of the 
electron-phonon coupling. On the other hand, in Ref. [2 ] decoherence has been derived even 
for zero temperature, where the bosonic phonon modes are in the vacuum state. 

Non-Markovian coupling of a system to a reservoir, has received increasing theoretical 
interest in recent years [22-31]. An important conceptual issue is to obtain a solution 
when decay rates attain negative values. In the literature, there exist numerical methods 
to take into account negative decay rates. These comprise quantum trajectory simulations 
in extended Hilbert spaces [32-36] and alternative quantum-trajectory methods [37, 38]. 
However, for presently studied physical systems, negative decay rates have been predicted 
only for extremely short time scales. 

We will show in this work that donor-based charge qubits, subject to the interaction with 
acoustic phonons, do naturally show this feature on time scales, that may be observable in 
experiment. Furthermore, the model developed here contradicts the result of Ref. [21], where 
decoherence is obtained even for zero temperature. From our work, the coherence time is 
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obtained as inversely proportional to temperature, leading to vanishing decoherence at zero 
temperature. 

2. Model of Phonon Scattering with a Localized Electron 

The Hamiltonian of the system under study is composed of three parts, 

H = H e +H v h + H e p h, (1) 

where the electronic part is in principle composed by the Bloch electrons plus the electrons 
bounded to the impurities. The quantized Fermionic field for electrons with spin projection g 
can be expanded as 

= £<WV«,q( r ) +E4,n,g<Mr), (2) 

n,q a,n 

where y„,q(r) and a ,n( r ) are the Bloch functions of itinerant electrons and the wavefunction 
of the electron localized at the donor site a, respectively. The operators c n ,<i,q and d a , n ,g 
annihilate an electron with spin projection g in the Bloch state l// n q and in the localized state 
flj „, respectively. In our case the impurities consist in pairs of donors (e.g. P) embedded in a 
semiconductor substrate (e.g. Si). The full electronic Hamiltonian reads 

n,q,g a,n,g 

+ \ H hA »-»' (d+,nd-,n' + dl in J+,n^ , (3) 
)i, n' 

where E n q and he a n are the energies of Bloch and localized electron, respectively. Here we 
added the third term which is responsible for tunneling between localized states at different 
donor sites with rates A n „/. 

The phonon Hamiltonian corresponds to the energy of vibrational excitations of the 
crystal lattice of the host semiconductor (Si). These excitations are bosonic in nature and 
we can write, 

k,cr 



where the bosonic operator b^ a annihilates a phonon with wavevector k, polarization o 
(o = 1 , 2, 3) and energy hcck^. We will consider only acoustic phonons, as the larger energies 
of optical phonons make them practically unpopulated at room temperature and below. 

2.1. Electron-phonon interaction 

The electrons couple to the lattice vibrations through the third part of the Hamiltonian, 

H e - ph = e J J 3 rp e (r) j£u„-VV(r-R„)}, (5) 

where the sum is performed over the lattice sites. Here e is the elementary electric charge and 
p e (r) is the electron density. Furthermore, V(r) is the Coulomb potential of a single ion of 
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the lattice with R„ and u„ being the equilibrium position and quantized displacement of the 
ions, respectively. 

In terms of the bosonic phonon operators we can write Eq. (5) as 

4-ph = U d 3 rp e (r) £ £ K k , K ,a (h,o + P_ k a ) ^ k+K > - r (6) 
y J k,cr K V 7 

where Y is the volume under consideration and the k and K sums are over the first Brillouin 
zone and the reciprocal lattice, respectively. The coupling constant is derived as 



Nh 

Kk.K.cr = iej- (k + K) -Ck^Vk+K- (V) 

2m £0k,(T 

with mo and N being the ionic mass in the unit cell and the number of unit cells in the 
semiconductor crystal. Here the Fourier coefficients of the ion potential are defined by 

y ( r ) = ^2>e' k - r . 

y k 

Using Eq. (2), the electron density can be written as 

p e (r) = p e (B) (r) +p e (BL) (r) +p e (L) (r), (8) 

where the different terms correspond to pure Bloch electrons (B), cross terms between Bloch 
and localized electrons (BL), and pure localized electrons (L). Inserting Eq. (8) into (6), 
the second term (BL) will produce phonon-induced transitions between Bloch and localized 
electrons. However, for T < 300K the thermal energy of phonons is not sufficient to drive 
such transitions. Therefore, we can neglect the BL term together with the B term, since under 
these circumstances the latter does not affect the dynamics of the localized electron under 
concern. Thus, we are left with the electron density 

Pe (L V) = EE^^ffflUWMr), (9) 

a,b n,n' 

For sufficient high temperatures, where the average thermal phonon energy, kgT ', is of the 
order of the electronic excitation energies of individual donor sites, phonon-induced electronic 
transitions n — > n' n') within the individual donor sites (a = b) or between the donor sites 
(a ^ b) may occur. In this case, resonant phonon scattering occurs that leads to a decoherence 
mechanism with a strong dependence on the inter-donor distance [39,40]. 

In this work we consider a different temperature range, i.e. kgT is much smaller than the 
electronic excitation energies of individual donor sites. In this case only the electronic ground 
states of the donor sites can be populated and the density (9) is further simplified to 

p e (L) (r)=£j»;(r)^(r), (10) 

a.b 

where d a = d a ,o,g with g being the invariant spin projection of the single electron and 
a (r) = (f) a ,o(r). Inserting Eq. (10) into the interaction Hamiltonian (6) we obtain 

#e-ph=^ £ ^K k , K ^F c , b (k + K)dy b (b k , a +b\ a ), (11) 
y a,b=±k,o K 
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where the form factor is defined as 



F ajb (k) = J ^r0*(r)^(ry*- r . (12) 

Given that the rms spread of the localized electron wavefunction is much larger than the 
periodicity of the semiconductor lattice, i.e. it spreads over several lattice sites, the form factor 
(12) is non- vanishing only within the Brillouin zone, F a &(k + K) « if K ^ 0. Furthermore, 
given that the wavefunctions of the two donor sites do not overlap, F a ^(k) = for a ^ b. 
Considering the medium as being approximately isotropic the interaction Hamiltonian (11) 
can then be simplified to 

#e-ph = ^Y,L K kFa(k)4d a (b k + P_ k ) , (13) 
T a k V 7 

where F a (k) = F aa (k) and where b k are annihilation operators of acoustic phonons with 
longitudinal polarization. The simplified coupling constant is 



Nh 

"^'W 4 - <14) 

with linear dispersion 0)^ = sk, sound velocity s, and with V k being the Fourier components of 
the ion potential. 

2.2. Electronic pseudo-spin 

In the basis of the localized-electron states \a) {a = ±) with (r\a) — a ( r ) an d f° r on ly a single 
localized electron being present, i.e. dld a \b) = 8 a ,b\b), one may define the electronic pseudo 
spin- j operator, 

&=^J»(&|, (15) 

1 a=£b 

S Z =~£«|a)H ( 17 ) 
^ a 

and the electronic identity operator 
4 = £|«)R 

a 

With these definitions the sum over donor sites in Eq. (13) is rewritten as 

Y,F a (k)dld a = S z f k + I e F k , (18) 

with 

F k = ^F a (k) 1 (19) 
/ k =F+(k)-F_(k). (20) 
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The term proportional to the identity operator / e can be shown to lead to a phonon- 
induced shift of the electronic ground-state energies he a Q of the two donor sites (a = ±), see 
Appendix A. It can be accounted for by defining the displaced phonon operator 

«k = £k + /3k, (21) 
with the displacement being given in Eq. (A. 3) and by redefinition of the electronic energies, 
see Appendix A. In this way the electronic Hamiltonian (3) can be simplified to 

6 e = h(OoS z + hAS x , (22) 

where ojq is the transition frequency between the donor-site electronic ground states, modified 
by the phonon-induced energy shift, cf. Eq. (A. 6) in Appendix A, and A = Ao,o is the tunneling 
rate between the ground states of the donor sites. The interaction Hamiltonian (13), on the 
other hand, becomes 

#e- P h = (gkd k + gkflk) , (23) 

k V ' 

where we defined the interaction rate 

K kf-k 

(24) 

2.3. Interaction rate 

Assuming a Yukawa potential for the screened ions of charge Ze, Eq. (14) becomes 



Nh k 

" l £ V2moffl|fctf+*?' ( 
with k s being the inverse screening length of the ionic potential. For semiconductors (Si) and 
T < 300K the wavelengths of populated phonon modes are much larger than the screening 
length, i.e. k <C k s , so that we may approximate 



k ~i Ze2 Nh( ° k (26) 
k Z 5£o^y 2m ' 

The ground-state wavefunctions at the two donor sites are displaced by the inter-donor 
distance vector d, i.e. 



1 

0±(r) = ; exp 

7CR 3 ± 



|r=Fd/2| 
R± 



(27) 



where R a (a = ±) are the Bohr radii of the s-wave orbitals of the two donor sites. Thus the 
form factors become 

±ik-d/2 
[l + (kR±/2)< 

so that together with Eq. (26) the interaction rate (24) becomes 



F ±( k ) = n ■ n. n 7^m , (28) 



D I hco k f e-™l 2 J™' 2 

8k hs y 2 Pm r \[i + (kR + /2) 2 } 2 [i + (*r_/2)2]2 ; 

where p m and "Vq are the mass density and unit-cell volume, respectively, and the deformation 
constant results as D = Ze 2 / (eQ^k 2 ). 
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3. Reduced Electronic Dynamics 



We assume here that during the free evolution of the system, tunneling is inhibited by either 
an applied potential barrier between the donor sites, or by the different electronic ground- state 
energies of the donor sites, |coo| 3> |A|, which can be provided for by the application of a DC 
electric field. Under these circumstances the Hamiltonian is given by 

k k 

3.1. Electron-phonon eigenstates 

We look for the eigenstates of the Hamiltonian (30), 

H\E)=E\E), (31) 
where from the Hamiltonian it is clear that they have the form 

\E)->\E m ) = \m)®\& m ), (32) 
with S z \m) =m\m) and with |<J> m ) being a solution of the phononic eigenvalue problem 

\&m), 03) 

with the spin-projected phononic Hamiltonian being 



H m = ^hcOkd^d^ + hm 

k 



~k 

This Hamiltonian can be diagonalized as 



£0D + L(^k + ^k«k) 



(34) 



H m = Y< hG)kdL \t.m^,m + hmCOQ, (35) 
k 

where we used m 2 = 1/4 and discarded a constant energy term. Here the operators 

«k,m = «k + ^0£k, (36) 

are phonon operators with displacements conditioned on the electronic state m. Their 
displacement amplitudes are given by 

«k = — ■ (37) 

(Ok 

The displacement of the phonon operators can be realized by the single-mode displacement 
operators 

D k (a) = exp (aal - a*d k ^ , (38) 

as 

dk, m = «k + ma k = E>1 (mat) d k D k (ma k ) . (39) 
Thus, the Hamiltonian can be written as 

H m = D f {{ma k }) H m D {{ma k }) , (40) 
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where 

H m — £ hco^d^dk + hmCQo, (41) 

k 

and where the multi-mode displacement operator is defined as 



D({m(Xk}) = Y\D k (ma k ) = exp 

k 



£m ^CCk«£ — OC^Ok 



(42) 



k 

As a consequence, the phononic eigenvalue problem can be rewritten as 

H m \® m )=E m \® m ), (43) 

where the transformed eigenstates are 

|<J> m >=D({ma k })|<J>,„>. (44) 

Given that H m is the Hamiltonian of an infinite set of harmonic oscillators, these eigenstates 
are identified as those of the harmonic oscillators, i.e. 

|$ m )-H{Ak», (45) 

with N k being the number of phonons in modes k. Consequently the energy spectrum is 
discrete and reads 

E m ->■ E m {Nk} = Y,nco k N k + hmcoo. (46) 

k 

Transforming back to the original frame, the complete eigenstates are therefore given by 
displaced number states: 

\Em,m) = \m) ®& ({ma k }) \{N k }). (47) 

The general solution of the reduced density operator of the spin system, i.e. traced over 
the phonons, reads therefore 



Ps(0= £ \m){m'\ £ £ p, 



m,{N x };m' : {N' x } 



({N k }\D({m'a k })tf({ma k })\{N k }) 



x exp 



it / 

J [ E {N k },m-E{ N £, m > 



(48) 



where we used for the phononic trace Tr(|y/)(0|) = (0|y/) and where P m) {N^};m',{N'} 1S me 
initial density matrix of the complete spin-phonon system. It can be easily seen, that the 
diagonal matrix elements of the solution (48) are constant, {m\ps(t)\m) = (m\ps(0)\m). 

We note that the initial reduced electronic density matrix elements result from Eq. (48) 

as 

(m\p s (0)\m) = £ p m ,{;v k }; m ,{jv k }, (49) 
(g\ps(0)\e) = £ £ PgWMKyfiKhm < 50) 
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and (e|/$s(0)|g) = (g|Ps(0)|e)*. The non-diagonal reduced matrix element (50) contains the 
Franck-Condon transition amplitude [41-43], well known from molecular physics, 

f{KUN*} = (iK}\D({oc k })\{N k }). (51) 

The presence of this factor is due to the fact, that the initial density matrix of the complete 
spin-phonon system, p m .{Ar k } ;m '.{A^}, is defined in the basis \E m ^ Nk j), that differs from the 
basis \m) ® | {A^k}) by the displacements of the phonon state. 

3.2. Preparation of the initial electronic state 

The preparation of a coherent superposition of electronic states \g) and \e) must start from a 
well defined state, which we assume to be the ground state \g). This can be experimentally 
realized by applying a voltage U ^> k^T between the donor sites. When thermal equilibrium 
is reached after the switching on of the voltage, the occupation probability of state \g) is 
very close to one. The coherent electronic superposition can then be achieved in two ways: 
Either by switching off the voltage to allow for coherent tunneling between the donor sites, or 
by applying THz electromagnetic radiation to perform a coherent Raman transition between 
states \g) and \e) [44-46]. To observe the decoherence of the prepared superposition state 
without the detrimental effects of coherent transitions between the donor sites, the coherent 
tunneling must be suppressed. This may be achieved either by applying an electric field 
creating a potential barrier between the donor sites or by switching on a sufficiently large 
voltage between the donor sites. 

As outlined above, the initial state p (0) of the complete system is the result of a state- 
preparation process in which the spin originally starts in its ground state \g) and is coherently 
transferred into a superposition of ground and excited states. In the preparation process, a 
fraction | ^ | 2 of the population is coherently transferred from the ground state \g) to the excited 
state \e). However, as in the excited electronic state each phonon mode receives a coherent 
displacement a^, cf. the eigenstates in Eq. (47), Franck-Condon type transitions [41^-3] 
occur for the phonons. 

The probability amplitude ^{Mk} to en d in state \E e t My \) is given by 

¥e,{M k } = £ £ ( E e,{M k }\S+\E gj{Nk} )(j)^ {Nk} 

= € E /{M k },{iv k }0g,{iv k }> (52) 

where ^{a^} are the probability amplitudes for starting originally from states \E g r Nk \). On 
the other hand the probability amplitude Vg^A^} t° stay in state \E g ^ Nk j) is 

Vg,{;v k } = V^-IWftW- (53) 
Thus the prepared complete initial state can be written as 

I V(0) > = £ {Yg,{N k }\E g , m ) + Ve,{N k }\E e ,{N k })) , (54) 
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which by insertion of Eqs (52) and (53) becomes 



|V(0)>= L <i>g,{N k } 



X 



1 - l^l 2 l^g.{iv k }) +^D({a k })\E ej{Nk} ) 
The corresponding initial density operator is p (0) = | i/a(0)) ( i/a(0) | and reads 



(55) 



p(o) =EE ^.wfe) (56) 

WW 



(i-l^l 2 ) l^.w^wl 
+ |^| 2 D({a k })| J E e ^ k} )( J E,. W} |D t ({a k }) 



+ rV 1 -|«l 2 K{At}><^«l fit ({ fl k}) 



+^V 1 - l^| 2jD ({«k})|£ e ,{yv k })(%<} 

We assume that originally the spin-phonon system started from the spin ground state and 
the phonons being in a thermal state, i.e. 

P = E P {N k }\ E 8,{N k })(E g ,{N k }l (57) 
W 

with 



P W 



Z-'exp ( -£ftiV k ) , Z= £ exp ( - J>iV k ] , (58) 
where fa = ftofy/ (k^T). Thus we have to substitute 

^, wc w ~* p w n 3 ^k ' (59) 

k 

by which the initial prepared spin-phonon density operator (56) becomes 

0(0) = E (6°) 

W 

x {(l-|^| 2 )|£ g ^ k} )(^ k} | 

+ |^| 2 D({a k })| J E e ^ k} )( J E e ^ k} |D t ({a k }) 

+ §* a/ 1 - \S\ 2 \E g , {Nk }){E e , { N k }\DH{a k }) 



, W I 

In the basis of the states (47) the corresponding density matrix reads then 

PftWaW = " l?l 2 )%}II^ k ' (61) 

k 

Pe,{N k };e,{Nl} = |€l E / '{M k }/{A' k },{M k }/{V k },{M k } ' ( 62 ) 
{^k} 



Pg.W^W = 5 V 1 " \^\ 2p {N k }f{ N '},{N k }- (63) 



Analytical model of non-Markovian decoherence in donor-based charge quantum bits 1 1 

Furthermore, the reduced initial prepared spin density operator becomes 
P5(0) = (l-|^| 2 )|g)(g| + |^| 2 |e)(e| | 

(64) 



(65) 



+ ^yi-|^|g}(,| + ^ v /l-|^|e>(g|, 
which corresponds to the perfectly pure superposition state 



Thus the presence of the phonons does not prevent a coherent preparation of the initial spin 
state. 

3.3. Spin Decoherence 

From the general solution (48), the time-dependent off-diagonal matrix elements read 



(g\Ps(t)\e) = e iaot ^ E Pg,mMK} 

mm 

x ({<}|D({a k })|{iV k })exp 
Inserting the initial state (63) this matrix element becomes 



-iY,o) k (N k -Nl)t 



( g \ps(t)\e) = e^'^Jl-\^Z 



12-7-1 



x Tr 



D t ({a k })D({a k (0})exp -£j3A 



where 



D({a k (t)}) =exp iY,G) k N k t D({a k })exp 



with 



-ia k t 



a k [t) = a k e 
The displacements can be joined to obtain 

ff({a k })D({a k (t)}) = D({a k (t) - a k }) 



x exp 



-iV |a k | sin(c<v) 

k 



so that Eq. (67) becomes 



(g\ps(t)\e) = e^'^^l-\^-Z- 1 

x Tr D({a k (t) -a k })exp ( -£j6A 

x exp — i^ \a k \ s'm(cQkt) 



(66) 



(67) 



(68) 



(69) 



(70) 



(71) 
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The thermal averages can be calculated in phase space as 



Tr 



D(a)e-P N 
Thus Eq. (50) reads 

where the time-dependent phase reads 
0(?) = £|a k |sin(fi)^), 

k 

and the decoherence rate is given by 

7(0 = 2£^£^sin(M 



j d 2^-a^ e -m\ 2 



I«| 2 /j8 



Sd^e-m 2 

i(Oot-i^(t) - f dxy(x) 
Jo 



k 

2k B T 
h 



£|oc k | 2 sin((0^)- 



(72) 



(73) 



(74) 



3.4. Decoherence rate 

The mode sum is evaluated and after integration over the spherical angles of the mode's 
wavevector, the rate becomes 



yit) 



k B TD 2 



dk\^ 



ks'm(skt) 



+ 



(2;r) Wp m i-oo [i + (kR a /2) 2 } 4 

cos(\d + st\k) — cos(|J — st\k) 



d[\ + (kR_/2) 2 ] z [\ + (kR + /2) 2 
It can be written as 

r(t)= EM0> 

a^b 

where we defined 

3 



MO = S 

'la 



+ 



1 



3 \T d r} a 

2 \ 2 

il 



i 

+ - 



i 

+ - 



il - il 



la 



ll - ll 



1+1 



1 — 



2 V^W 4 



+ 



+ 



la ll-$ll ~ 

2 ll-ll. 



-2t/{x d r\ a ) 



2\l+t/X d \/T} a 



-2\\-t/T d \/r\ a 



(75) 



(76) 



(77) 



where we defined the relative Bohr radii 
T) a =R a /d, (a = ±), 

and the time needed for a phonon to travel between the donor sites, x d =d/s. The temperature 
dependent rate is defined by 

T 
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where (O c i = 2k/ is the angular frequency of a phonon with wavelength X =d. The material 
dependent temperature scale is defined for convenience as 



where Nj is the number of unit cells within the volume d . 
4. Discussion 

4.1. Time-dependent decoherence rate 

The initial value of the decoherence rate can be shown to be y(t = 0) = 0. Moreover, it 
is easily shown that it asymptotically decays to zero for large times, y(t — > °°) = 0. The 
dimensionless and temperature-independent decoherence rate y(t) /Tj is plotted in Fig. 1 as a 
function of the dimensionless time t/tj. Here the mean relative Bohr radius r\ = Y* a la/2 has 
been chosen as rj =0.1, whereas the relative difference of (relative) Bohr radii, o = Arf/r] 
with At] = r] + — T] , is varied for the three different curves. Curves for relative differences 
of Bohr radii smaller than 0. 1 can be shown to converge to the blue curve with a = 0. 1 . The 
positive and negative peaks become more pronounced for decreasing ratio of Bohr radius over 
inter-site distance, see Fig. 2. It is observed that whereas the widths of the peaks decrease, 
their heights increase. 

It is readily recognized that the maximum decoherence rate always occurs at the time 
t ~ T d mm(r} a ), which is the time required for a phonon to cross the smallest donor site. A 
negative peak of the decoherence rate, i.e. a maximal recoherence rate, is observed always 
around the same time t ~ %d- This corresponds to the time needed for a phonon to travel 
between donor sites. Thus, the timescales can be interpreted as the time needed for the 
phonon to either travel within an impurity site or between the sites. In the first case, when 
the phonon (as quasiparticle) is traveling within one of the sites it destroys the inter- site 
coherence and does this at maximum rate at the moment it has propagated a distance equal 
to the smallest Bohr radius of the impurity states. In the latter case the negative minimum, 
i.e. the recoherence peak, appears when the phonon has traveled the distance between the 
impurities, thereby correlating their states. 

The first time scale of the positive peak is produced by the first term in Eq. (77), whereas 
the second timescale of the negative peak comes from the third term. The second term, on the 
other hand, compensates the third one for times t <C x c \ and provides for y a ^{t = 0) = 0. It 
can be easily seen that when expanding the decoherence rate (76) with respect to the relative 
difference of (relative) Bohr radii o we obtain 



Here yb(t) is the decoherence rate for the limiting case of equal Bohr radii, o = 0, which is 
obtained as 




(78) 



y(t) = y (t) + G (a 2 ) . 
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Figure 1. Dimensionless decoherence rate y(/)/rY as a function of the dimensionless time 
t/Xd for tj = 0.1 and relative differences of the (relative) Bohr radii being: a = 1.0 (red 
curve), a = 0.5 (green curve), and a = 0.1 (blue curve). Curves with a < 0.1 become 
indistinguishable from the blue curve. 




Figure 2. Dimensionless decoherence rate y(f)/rY as a function of the dimensionless time 
tjx t i for a = 0.1 and mean relative Bohr radius being: TJ = 0.1 (red curve), rj = 0.08 (green 
curve), and rj = 0.05 (blue curve). 
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The fact that the decoherence rate may attain negative values clearly indicates a strong 
non-Markovian behavior. Thus, a certain memory can be assigned to the phonon reservoir, 
with its timescale being determined by the time needed for a phonon to travel the distance 
between the impurities. The occurrence of negative decoherence rate and non-Markovian 
behavior in an exactly solvable system, such as presented here, is of high interest from 
the viewpoint of quantum-trajectory theory. There, a controversy exists on how to find the 
correct trajectories when decay rates may become negative [33,36-38]. However, up to now 
physically relevant models where such peculiar rates can be observed at relevant timescales 
are largely unknown. Thus, our model may be the first to show such features as experimentally 
observable effect. 



4.2. Dependence on temperature 



The electronic coherence, i.e. the off-diagonal electronic density matrix elements, decay with 
the function 



exp 



- fdt'y{t') 
Jo 



cf. Eq. (73). It can be rewritten as 

g(t) = [G (t/x d )] T/T \ 
The decay function at T = Tq with dimensionless time x = tjx^ is obtained as 

2% rx 



Gq(x) = exp 



Tr Jo 



dxy(x d x) 



(79) 



(80) 



Performing the time integral for the case o = 0, i.e. for identical donor sites, one obtains 



lnGoOc) = — —A 

+ 0(1 -jc) 



+ B 



B 



1 +x 



1 
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+ [1-0(1-*)] 



B 



B 
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+ B 



25(0) 



(81) 



where 0(*) is the Heaviside step function, whereas A(x) and B(x) describe the effects of the 
positive and negative peak of the decoherence rate, respectively, 
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Figure 3. Temperature dependence of the electronic decoherence function g(t) as a function 
of dimensionless time t/Tj for 77 =0.1, <7 = 0.0 and temperatures: T /To = 0.1 (red curve), 
T /Tq = 0.02 (green curve), T /To = 0.01 (blue curve), and T /To = 0.002 (purple curve). 

Note, that the temperature-independent function (80) depends in general on 77 and a, but it 
does not depend on any parameter of the semiconductor material, such as sound speed etc. 
This feature is provided for by the proper definition of the critical temperature in Eq. (78) and 
the use of t/Xd as dimensionless time in Eq. (79). 

In order to get a grasp on the experimental parameters involved, we consider the donors to 
be phosphorus atoms, separated by a typical distance of the order of d ~ lOnm, embedded in 
a silicon substrate. Silicon has a mass density p m = 2.33 x 10 3 kgm -3 , a deformation constant 
D = 8.6eV, and sound speed s = 9 x 10 ms . With these values the time for a phonon 
to travel between the donor sites is of the order of magnitude ~ lps. Most importantly, 
the characteristic temperature scale, as defined in Eq. (78), is of the order of Tq ~ 10 3 K for 
this material. In consequence, all relevant operational temperatures are much lower than this 
temperature scale (T <C Tq). 

In Fig. 3 the decay function (79) is shown for r\ = 0. 1 (a = 0) and various values of T 
with T <^Tq. It can be readily seen that apart from deviating from the Markovian exponential 
decay it additionally is a non-monotonic function, i.e. it first decays and then rises again 
due to the negative decoherence rate at t ~ Xd- A simple definition of a coherence time x c 
by requiring g(x c ) = l/e, being appropriate for an exponential decay, leads to problems of 
multi-valuedness in this case. 

It is more convenient to specify the transient mean decoherence rate 7, corresponding to 

the rate averaged up to the time Xd when the recoherence starts, 

1 r T d 
f=— dty{t). 
Jo 
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Employing the results in Eqs (79) and (80), and using r}T,j = R/s, the order of magnitude of 
the coherence time of the system, x c = 1/ y, is obtained for o = and r\ <C 1 as 



Thus for T = OK the coherence time becomes infinite and no decoherence occurs^. This result 
is in contrast with recent results using a Wigner-Weisskopf calculation but truncating the set 
of coupled differential equations, where decoherence can be observed even for T = OK [21]. 

5. Summary and Outlook 

In summary, we have studied the decoherence dynamics of a charge donor quantum bit and 
have obtained results in agreement with numerical studies [20] but in disagreement with 
the zero-temperature theory proposed in Ref. [21]. Our treatment is entirely analytic and 
exact and provides an analytical framework of a strongly non-Markovian physical system 
for testing recently published quantum trajectory methods for systems with negative decay 
rates [25,33,36-38]. 

Our theory reveals that the dynamics of the electron-phonon system is dominated by 
Franck-Condon type transitions that entangle electron and phonon states by an electron 
conditioned coherent displacement of the phonon modes. As a consequence, the reduced 
electron dynamics shows two strong non-Markovian features. Firstly the initial decay is not 
exponential but quadratic and secondly a later recoherence is observed due to the occurrence 
of a negative decay rate. This recoherence leads to a stationary and appreciable value of 
electronic coherence for temperatures T <C Tq, with Tq ~ 1000K, which thus should be 
observable at room temperature. Furthermore, as the dynamics of the electronic coherence 
is described by an analytic function, its dependence on temperature is obtained in analytic 
form as a simple scaling law. 

In conclusion, our results indicate that charge donor-based implementations of quantum 
bits are promising candidates for quantum memory and computation at feasible temperatures. 
However, further investigations are needed to study other types of decoherence mechanisms 
that may affect this type of system, such as coupling to nuclear spins or electromagnetic 
fluctuations. Furthermore, the effects of coherent tunneling could also be included in the 
future in a perturbative manner. 
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Appendix A. Phonon-induced shift of electronic energies 

After inserting Eq. (18) into the Hamiltonian (13) and using Eq. (4), one obtains the complete 
Hamiltonian as 

H = He + Y,nca k bib k + ^Y< K k{Szf(k)+N e F(k)) (S k + Sl k ) . (A.l) 

k k 

Given only one localized electron, N e — > / with / being the (electronic) identity operator, so 
that the phonon part can be diagonalized by the use of the shifted phonon operator 

d k = b k + p k , (A.2) 

with the shift being 

Performing the approximations as described in Sec. 2 the electronic Hamiltonian reduces 

to 

H e = h(£ +7 o-£-,o)S z . (A.4) 

Due to the use of the displaced phonon operator (A.2), this electronic Hamiltonian is modified 
to obtain the complete Hamiltonian as 

where the shifted electronic transition frequency is 

1 . 

W 



ojo = (e+,o - e-,o) - ^ £ **/(k) {Pk + P *k) • (A.6) 
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